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ABSTRACT 


The space and time continuity of numerical fronts as operationally 
produced by the Fleet Numerical Weather Facility, Monterey, California 
(FNWF), are investigated at 1000, 850, 700, 500, and 300 mbs for the 
four-day period January 25-28, 1967. Front locations, intensities, 
movements and slopes are examined over North America and compared with 
those from other analysis centers. Vertical structures of baroclinic 
zones and their changes with time are also investigated. Finally, a 
study is made of the FNWF computer-processed temperature fields compared 
to manually analyzed temperatures. 

Results indicate reasonable vertical consistency of frontal infor- 
mation in the central and eastern United States below the 700-mb level. 
Time continuity is best maintained at the 700-mb surface. A fictitious 
displacement of numerical fronts toward the warm air is observed at each 
of the mandatory pressure levels considered. The temperature fields of 
FNWF are found to produce smoothing of thermal perturbations associated 
with open frontal waves while relaxing the thermal gradient which leads 


to widening and weakening of baroclinic zones. 
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A CASE STUDY OF THE SPACE AND TIME CONTINUITY 


OF NUMERICAL FRONTAL ANALYSIS 


l. Introduction 

The need for a numerical method of frontal analysis has been recog- 
nized for some time. In the search for a front-location parameter to 
which numerical techniques could be applied, Renard and Clarke (11) chose 
potential virtual temperature (0). Specifically, the parameter is de- 
fined as the directional derivative of the gradient of potential virtual 
temperature along its gradient on a constant pressure surface, hereafter 
symbolized as GGO. Thus, 

° = 

vive! vo = M VIVO| q No 


GGO=- 
(ve (1) 


where No iS a unit vector in the direction of VO. 

Figure 1 shows the relationship of 0, VO, and GGO. Areas between 
lines of GGO maxima (i.e. ridges) and minima (i.e. troughs) represent 
hyperbaroclinic zones. The ridge line depicts the warm air boundary of 
the baroclinic zone and is thus the numerical analog of manually ana- 
lyzed fronts. 

The Navy Fleet Numerical Weather Facility, Monterey, California 
(FNWF), operationally computes a hemispheric field of. GGO on a twice- 
daily basis. The product is available to major weather centrals and 
facilities as an aid in frontal analysis. The Central Analysis Office, 
Meteorological Branch, Department of Transport, Montreal, Canada (CAO), 
is experimenting with a similar approach to objective frontal analysis; 


however, this thesis study is primarily concerned with the FNWF method. 


ISuch zones will be referred to as baroclinic zones in subsequent 
paragraphs. 
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Since 9 is derived in a unique way, it is of interest at this point 
to summarize :briefty the FNWF automatic processing technique. Surface 
observations and mandatory-level radiosonde data plus extrapolations from 
ship reports and selected aircraft observations comprise the raw data 
field. Using these data a static stability parameter is derived for five 
layers (1000-775 mbs; 775-600 mbs; 600-450 mbs; 450-300 mbs; and 300-200 
mbs). The stability parameter is a constant for each of these layers at 
each radiosonde station, but is allowed to vary in the horizontal. 
Using the 1000-mb height, 1000-500-mb thickness and static stability as 
specified at intersections of a 63x63 linear grid, the 0 s at mandatory 
levels may be computed. The Q's so determined (hereafter called pro- 
cessed 9) are the basis for calculation of GG@ in equation 1:. De- 
tails of this processing technique may be found in (1). 

It is the purpose of this thesis to investigate the space-time 
continuity of baroclinic zones as indicated by the field of GGO and re- 
lated analyses for a winter synoptic situation, 25-28 January 1967. 

The ultimate goal is to enhance the usefulness of the GG@ analyses and 
prognoses as received by the Naval Weather Service operational consumer. 

The research comprises three phases of investigation and attempts 
to answer the specific questions posed below: 

a. Vertical Consistency 

(1) What is the geometric relationship of the tropospheric 
baroclinic zones at each of several mandatory levels? 

(2) How does the orientation and slope of numerical fronts 
compare to their manua lly? drawn counterparts? 
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front. 


"Manually drawn is to be considered synonomous with ''conventional'' 


(3) Is there vertical consistency in position and intensity of 
GGO centers as found along the GG0 ridges? 

b. Time Continuity 

(1) What is the relationship of the baroclinic zone from one 
synoptic time to another at representative mandatory levels and how do 
baroclinic layers change with time? 

(2) Do the position and intensity of 669 centers exhibit time 
continuity? 

(3) Is the movement of numerical fronts consistent with the 
horizontal wind field and movement of manually analyzed fronts? 

(4) Is there a usable relationship between movement of pressure 
centers and GG0 centers? 

(5) Is there diurnal variation of movement and intensity of 
GG0 centers? 

c. Relative Comparison of Various Frontal Analyses 

(1) How do FNWF frontal positions compare with CAO and the 
Environmental Sciences Services Administration (ESSA) locations? 

(2) What are the existing and potential merits and deficiencies 
of the FNWF frontal analyses, the knowledge of which would aid the con- 
sumer in the operational use of GGO charts? 

As an adjunct to the primary objectives, a kinematic geostrophic 
expression for computing frontogenesis (10) is applied to the O field. 
Employing the geostrophic wind, the working equation may be symbolized as: 

E 
FG=-V(%- YS )-Ve=-v(%-7,)-Ve 


I vel (2) 


Positive (negative) values indicate an instantaneous increase (decrease) 


in horizontal gradients of 9. 


2. Area, Data, and Charts 

The region studied for vertical consistency and time continuity was 
limited to the Northern Hemisphere south of 55% latitude and between 652 
and 125º west longitude. For that portion of the investigation dealing 
with comparison of frontal positions the area was extended as indicated 
in Section 5. 

Computer printouts of numerically contoured 0, Vo, GGO and FG at 
standard tropospheric levels (1000, 850, 700, 500, and 300 mbs) for the 
synoptic times 1200 GCT 25 January to 0000 GCT 28 January 1967 were ob- 
tained from FNWF. Each printout is a 63 x 63 grid-point field on a 
polar stereographic projection with a scale of 1:30,000,000. The mesh 
length is 381 kms. GGO ridges and troughs were drawn by hand. Any 
positive GGO zone with less than six grid-point values equal to or 
greater than +0.05 C/(100 km)? or which did not contain one value of at 
least +0.15 C/(100 km)? was ignored as meteorologically insignificant. 
ESSA's Northern Hemisphere surface maps and upper air analyses for the 
United States and Canada were used for comparison to FNWF charts. 

During the period of study an intense storm developed from a 
frontal perturbation located in the Texas Panhandle area 25-26 January 
1967 and deepened as it moved northeastward through the central states 
(Figures 2-7). The situation was a classical example of a cold air 
injection (5) over Texas with subsequent deepening of the primary low. 
While hand analysis of frontal positions must of necessity be used as a 
basis of comparison, it was believed, apriori, that the situation east 
of the Rocky Mountains would show a minimum of disagreement among ex- 
perienced meteorologists in locating fronts (see Figures 8-10). Frontal 


positions over the mountainous western United States are not nearly so 


well defined and continuity of conventional frontal analyses was more 
difficult to maintain. For this reason the focus of attention is on 
the central and eastern region and it is this area which generated most 
of the discussion and conclusions that follow. 
3. Vertical Consistency 
Thirty vertical cross sections of North American baroclinic-zone 
boundaries were plotted, twenty of which are shown in Figures 11-15. 
Vertical slices were selected in order to examine the baroclinic zone 
at approximately 5% of latitude either side of the peak of the warm 
sector (as shown by the ESSA analyses). Once these positions were deter- 
mined, additional sections were taken at approximately 10? intervals. 
This interval was shortened in a few cases in order to aid in evaluation 
of continuity. When fronts of short horizontal distance were encountered, 
a position near the mid point was selected. Each vertical slice was 
taken in a direction perpendicular, or nearly so, to the 850-mb GGO ridge. 
The selected positions are shown on Figures 16-21 by solid straight lines. 
With reference to the cross-sections of baroclinic-zone boundaries 
(solid curves in Figures 11-15), the cold (warm) side boundary is to the 
left (right) in each case. The width of the baroclinic zone at each 
pressure level is indicated in degrees of latitude at a position between 
the two profiles at the particular level. Similarly, the 0 values at 
each level considered are displayed in degrees absolute at the outer 
edges of the thermal zone. Values between levels represent frontal 
slopes coded as horizontal displacement of the frontal position (9 lat.) 
per 5000 feet change in height. Positive (negative) values represent 
slopes toward the cold (warm) air. Thus, 1.5 indicates a slope of ap- 


proximately 1 mile/105 miles. The dashed curves are vertical plots of 


ESSA fronts. The ESSA surface position is plotted at the 1000-mb level. 
ESSA analyzes fronts at 850 mbs but not at 700 mbs, therefore, tempera- 
ture profiles along the intersection of the vertical section and the 
pressure surface were plotted. The point of maximum curvature of the 
profile was used as the 700-mb frontal position. Sections have been 
arranged in groups to facilitate inspection of changes in time, but for 
the present, consider each section individually at a particular time. 

In general, the GGO trough profiles parallel the GGO ridges. The 
average width of the baroclinic zone (º latitude) at each level, for all 
thirty sections, is given below: 


1000 mbs 850 mbs 700 mbs 500 mbs 300 mbs 
6.3 6.8 6.6 Jak 6.5 


Potential virtual temperature is observed to increase with height. 
Average increases (in degrees absolute) through each layer are shown 
below: 


1000-850 mbs 850-700 mbs 700-500 mbs 500-300 mbs 
6.4 6.9 10.6 8.6 


Kreitzberg (8) states that broad baroclinic zones average 600 miles in 
width with 8° temperature change, while the narrower, intense zones 
average 130 miles in width with a contrast of 5°C horizontally across 
the zone. O'Connor (9) observed that in the cooler season the potential 
temperature of the polar front may increase from 2989K near the surface 
to 302°K at the 500-mbs level. The data presented above agrees reason- 
ably well with the findingsof Kreitzberg and O'Connor except that the 
widths of the majority of FNWF baroclinic zones appear to be biased 
toward the high values. The number of cases considered is insufficient 
to make any conclusions regarding changes in width with height, but the 


relatively narrow widths at 300 mbs were surprising. Inspection of the 
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tropopause heights during this period, however, indicates intersections 
with levels in the 500-300-mb layer and this makes interpretation of 
300-mb frontal information more complex. 

Geometrically below 700 mbs the baroclinic zone generally s lopes 
toward the cold air or is nearly vertical, while above 700 mbs both 
negative and positive slopes are numerous. Although negative slopes 
have been found at low levels, such occurrences in the free atmosphere 
are to be doubted (4). FNWF frontal slopes in 39 separate layers were 
compared with ESSA frontal slopes. Of these, 6 were identical; 18 were 
within 21.0°lat./5000 feet; and 21 differed by more than t1.0º1at./5000 
feet. 

Figures 16-21 are composites of GGO ridges at the standard levels 
for 1200 GCT 25 January to 0000 GCT 28 January 1967. Notice that over 
the continental United States eastward of central Texas the vertical 
consistency of numerical fronts is quite reasonable. Over the western 
United States, Mexico and Central America the slopes suggested are total- 
ly unrealistic. The frontal contours in the Gulf of Mexico region at 
1200 GCT 27 January and 0000 GCT 28 January demonstrate the numerical 
problems encountered as frontal systems move into an area of sparse data 
and/or when the thermal gradient weakens, while the western United States 
indicates well the fictitious frontal information introduced by extra- 
polation of data down to sea level. 

Centers of positive GGO (i.e. GGO maxima) at each level were plotted 
for each synoptic time. These plots are shown in Figures 22-27. In 
some instances GGO centers at 1000 mbs can be identified with maxima at 
higher levels, however, for the most part, there seems to be a poor 


relation. All centers observed at each level over the six synoptic 


times are combined in Figures 38-42 in order to study movement. A com- 
ment regarding the relative number of GGO centers at each level is per- 
tinent here. Total numbers observed at each level were 74, 48, 50, 40, 
and 61 at 1000, 850, 700, 500 and 300 mbs, respectively. The sharp de- 
crease between 1000 and 850 mbs is probably due to the effects of sur- 
face-induced heating and cooling and extrapolation problems in 1000-mb 
data. The increase at 300 mbs is more difficult to explain, but may be 
due to the fact that the tropopause complicates matters or that inten- 
sities of many relative maxima at this level are weak so that small 0 
errors would produce numerous centers. Should these hypotheses be correct, 
many individual maxima at 1000 mbs and 300 mbs have little synoptic scale 
meaning and correlations with adjacent levels are not to be expected. 

Intensity values of GGO maxima decreased with height in general. 
Greatest observed values were 7h, 61, 46, 49, and 52 (ºc/( 100 km)? 
at 1000, 850, 700, 500 and 300 mbs, respectively. 
4. Time Continuity 

The vertical cross-sections of Figures 11-15 have been arranged in 
five groups (A-E). Each plot in a group is a cross-section through a 
similar area of a frontal system at a different synoptic time. Thus, 
for example, by inspecting Group C an examination of time changes ina 
baroclinic zone at a position nearly fixed relative to an identifiable 
frontal feature can be made over a period of six synoptic times. 

Consider changes occurring in time through the layer 1000-700 mbs 
during the period 1200 GCT 25 January to 0000 GCT 28 January 1967: 

Group A - Sections in Group A were taken through a warm frontal 
region which became less intense with increasing time and moved from the 


eastern continental area to the western Atlantic Ocean after 1200 GCT 





25 January. The width of the baroclinic zone appears to increase from 
1200 GCT 25 January to 1200 GCT 26 January. During the period covered 
by Group A cross-sections, the nearest GGO maximum decreased from 65 to 
55 (0C/(100 km)? at 1000 mbs, increased from 30 to 47 at 850 mbs and 
showed no consistent pattern at 700 mbs. The numerical frontal slope in 
the layer 1000-700 mbs became steeper as did ESSA's frontal surface. 

Group B - This group is from a region analyzed by ESSA as a cold 
front until 1200 GCT 26 January and as a stationary front at 0000 GCT 
27 January. The baroclinic zone increased in intensity after 0000 GCT 
26 January. ESSA's frontal slopes remained near ly constant during the 
period while FNWF's became steeper until 0000 GCT 27 January. The near- 
est GG0 centers are the same as for Group A and were discussed above. 

Group C - Group C is a sequence of six different synoptic times. 
ESSA analyzed the frontal region as a cold front until 1200 GCT 26. January 
and as a warm front for the latter three synoptic times. No consistent 
pattern of change in width of the FNWF baroclinic zone or in the slope 
of ESSA or FNWF is noted. Adjacent GG0 centers fluctuated at 1000 mbs, 
increased steadily at 850 mbs until 0000 GCT 27 January, decreased at 
1200 GCT 27 January and then increased. The intensity fluctuated 
irregularly at 700 mbs. 

Group D - Sections in Group D were taken in an area of intense cold 
air advection. Winds in the cold air were on the order of 40 to 50 knots 
at 850 mbs and 50 to 70 knots at 700 mbs. The cold air injection was in 
progress at the time of the first section in the sequence. The thermal 
zone is narrow throughout the period. Slopes are observed to be steep. 


Group E - The frontal area of these sections was similar to Group D, 


but displaced to the west of the intense cold air injection. The geo- 


graphy of the area includes the mountains along the Texas-Mexico border 
and the western coast of the Gulf of Mexico. Irregularities in ESSA and 
FNWF frontal contours are to be noted. 

Above 700 mbs large space and time variations in the slopes of all 
baroclinic zones occur, possibly due to upper tropospheric fronts or jet 
streams. 

Again with reference to Figures 11-15, maximum variation of frontal 
Q values at specified levels, during the time intervals represented by 
each group, are shown in the following table (°C). It is to be noted 
that the least change occurred at 850 mbs and through a warm front sector 
(Group A). The former again highlights the utility of a free atmospheric 


level compared to the surface (or 1000 mbs) for frontal analysis. 


Level Group A Group B Group C Group D Group E 
1000 mbs 3 L 6 6 3 
850 mbs 2 2 10 3 2 
700 mbs 2 4 4 5 7 
500 mbs 5 y 4 4 L 
300 mbs 8 6 3 4 7 


Composites of the numerical frontal positions in time at the manda- 
tory levels 1000-300 mbs are included in Figures 28-37. Time continuity 
1S maintained best at 700 mbs. Frontal patterns can be identified from 
one synoptic time to the next and movements are reasonable. The picture 
at 850 mbs is much improved over that at 1000 mbs, but some unrealistic 
behavior is noted, as over the Gulf of Mexico from 0000 GCT to 1200 GCT 
27 January. Such movement may be due to lack of data, erroneous reports, 
or the sudden appearance of correct reports where such were previously 


absent. In any event, this example demonstrates the possibility of 
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unrealistic movement in sparse data areas. Deviations of numerical 
positions from hand drawn locations are given in Section 5, but since 
deviations did not increase with time, it is concluded that the average 
speed of FNWF fronts is consistent with that observed for conventional 
fronts. This is not necessarily true for individual 12-hour intervals 
as compared to ESSA positions. Movement evaluation based on the mean 
wind perpendicular to the front was attempted but the results were in- 
conclusive. Time continuity diminishes at 500 mbs relative to the lower 
levels and is poor at 300 mbs. 

The plots of GGQ maxima at each level considered (Figures 38-42) show 
that continuity is about equal at 850 and 700 mbs. In cases where GGO 
center positions exhibited continuity from one synoptic time to the next, 
but the hand analysis did not depict a front in the area, cyclone tracks 
were plotted. In most cases a similarity in movement of GGO and cyclone 
centers was observed E speeds varied. Deepening of cyclones ap- 
peared to be accompanied by increases in intensity of GGO centers but 
the number of cases observed was insufficient to make any conclusions. 
5. Relative Comparison of Various Frontal Analyses 

A comparison was made of FNWF frontal positions with those located 
numerically by CAO and manually by ESSA for the period 0000 GCT 26 
January to 1200 GCT 27 January. The CAO method (3) uses the GGO para- 
meter as well as others and processing procedures differ from FNWF. 

The FNWF 1000-mb GGO analysis was not used in this portion of the 
study,. but rather a frontal analysis derived from the mean temperature 
in the layer 1000-700 mbs (10/7 model). Since a number of fronts is 
produced by both the FNWF and CAO methods which have no counterparts in 


ESSA analyses or orientations were such that no meaningful comparisons 
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could be made in some cases, the following critera were established as 
basis for consideration: 

(a) Length of numerical front must be at least 10º of latitude. 

(b) At least one point on the CAO front must be within 59 of 
latitude of FNWF position. When comparing with ESSA analyses one point 
on either CAO or FNWF fronts within 5º of latitude of ESSA position was 
sufficient for both to be compared. 

(c) No comparison was made if the average orientation of a front 
was such as to make an angle of more than 45? with its counterpart. 

Deviations were measured at intervals of one grid length along and 
normal to the FNWF front for FNWF versus CAO, while ESSA's frontal 
position was similarly used as the base for the FNWF/CAO versus ESSA 
comparison. Deviations were measured to the nearest $ degree of lati- 
tude. Average deviations and average algebraic displacements toward 
cold or warm air in degrees of latitude are given in Tables 1-4, 

Deviations shown in Table 1 and 2 reflect a comparison at 850, 700, 
and 500 mbs only since the 1000 and 300-mb positions were not available 
from CAD. The area considered encompassed the west longitudes of the 
Northern Hemisphere. In Tables 3 and 4, FNWF's 10/7 model analyses were 
compared with ESSA's surface positions; here the area was the entire 
Northern Hemisphere. FNWF/CAO comparisons with ESSA were made only at 
850 mbs and the area was the continental regions of southern Canada, the 
United States, northern Mexico and the adjacent coastal waters. The 
cumulative average deviations and algebraic displacements are based on 
the average at each synoptic time weighted in proportion to the number 
of points measured. 


Table 1 indicates that FNWF-CAO deviations increase with height and 
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Table 2 shows that the average CAO position is displaced to the cold air 
side of its FNWF counterpart, the displacement increasing with height. 
Table 3 indicates that on the average, CAO positions are closer to the 
ESSA analyzed locations than are the FNWF fronts at 850 mbs. The 1000-mb 
and 850-mb deviations of FNWF from ESSA are equivalent, but Table 4 shows 
that at 850 mbs FNWF fronts have an average algebraic displacement rela- 
tive to ESSA fronts of about 2.3º latitude toward the warm ". The 
average algebraic displacement of the FNWF 1000-mb positions and the CAO 
850-mb fronts approximates the ESSA positions. 

Because of the results in Tables 1-h and since GGO is utilized ex- 
clusively by FNWF and inclusively by CAO for frontal analysis, an exami- 
nation of the FNWF processing method was undertaken to see if this fea- 
ture accounts for the difference noted. 

As stated earlier, the FNWF processing procedure involves the cal- 
culation of a stability parameter for specific layers which is constant 
through the layer. This is equivalent to obtaining an average lapse 
rate through the layer. Such a procedure produces deviations from the 
actual temperature at specified levels. The magnitude of the deviation 
depends on the variability in the sounding and the location of extremes 
(i.e. inversions) relative to the level where temperatures are desired. 
Such a conclusion is supported by statistics gathered by Schardt (12) 
who compared FNWF processed (virtual) temperatures to the observed dry- 
bulb temperatures at 1000, 850, 700, 500 and 300 mbs. His data showed 
that deviations in the range +1.0 to 2.9 and -1.0 to -2.9°C occurred at 
39% of the grid points at 1000 mbs, 53% at 850 mbs, 42% at 700 mbs, 25% 
at 500 mbs and 46% at 300 mbs. Further, the data showed that 16%, 8%, 


2.54, 84 and 16% of the grid points at 1000, 850, 700, 500, and 300 mbs, 
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respectively, had processed (virtual) temperatures colder than observed 
dry-bulb temperatures and hence are unquestionably in error. Such errors 
and others not isolated in this analysis could easily produce spurious 
baroclinic zones or cause artificial displacements of numerical fronts. 
The likelihood seems more probable in regions of weak thermal gradients. 

The arguments stated above do not explain the displacement of FNWF 
fronts toward the warm air side of hand drawn fronts using actual data 
on constant pressure surfaces. Figure 43 represents a simulation of the 
computer calculations of GGO assuming a correct temperature field. It 
also assumes a situation in which VO is parallel to Vivol so that the 
cosine of the angle between VO and Vol is one, when taking the dot 
product in the GGO formula (equation 1). Three temperature profiles are 
shown, as based on "grid point" data. The solid arrow indicates the 
actual point of maximum curvature (GGO maximum) and the dashed arrow 
shows the probable point at which the FNWF numerical method would indi- 
cate a GGO maximum. The upper profile represents the most intense front 
with decreasing intensities in the lower two. Note that the numerical 
computations displace the GGO maximum toward the warm air in all cases. 
Notice also that the numerical position would not be altered even though 
the dashed curve in the upper profile existed. Thus, an error in posi- 
tioning the front ranges from approximately % to I grid length (i.e. 1-3? 
latitude). The same feature is observed relative to the GGO troughs ex- 
cept the displacement is toward the cold air. The combined result could 
artificially broaden the baroclinic zone by as much as six degrees of 
latitude, 

In the two lower profiles where the change of temperature was less 


intense, the numerically indicated position approaches more closely the 


Q 
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actual position. The difference would be greater, however, if the point 
of maximum curvature occurred farther from the grid point. 

A further attempt to examine observed versus processed temperature 
fields was made by converting the ESSA analyzed dry-bulb temperature 
field to potential temperature and superimposing the FNWF processed field. 
Figures 44-46 show three such examples. To be sure, moisture effects have 
not been taken out of the FNWF field, but it is readily apparent that per- 
turbations of a synoptic scale (open frontal waves) are smoothed, the 
gradient is relaxed and regions of maximum change of gradient are dis- 
placed toward the warm air in the FNWF field. 

Because of the apparent displacement of FNWF fronts toward the warm 
air, a search was made for a correction which could be applied to the 
numerical frontal analyses at the mandatory levels. An arbitrary correc- 
tion toward the cold air of one-half the distance from the GG9 ridge to 
the 0267 (100 km)“ GG0 isoline was made and another comparison with CAO 
and ESSA was undertaken. Tables 5 and 6 show the results. Notice that 
the deviations of FNWF from CAO positions were reduced at every level; 
also the FNWF position relative to ESSA locations improved. Applying the 
adjustment to the FNWF 10/7 model positions resulted in greater devia- 
tions than were observed without the correction and consequently are not 
shown here. 

The 0.0 C/(100 kms ) 4 isoline is not computer plotted on FNWF GGO 
charts, but the +0.05 C/(100 kms)Žisoline is plotted; thus a displacement 
of the GGO ridge one-half the distance to the latter isoline will accom- 
plish near ly the same result. However, such a correction should not be 


made to the frontal product generated by the 10/7 model. 
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6. Frontogenetical Properties 

The applicability of equation (2) as an indicator of frontogenesis 
(FG) was investigated at the 850-mb and 700-mb levels in an area roughly 
comprising the continental United States. The average values of FG at 
each grid point for two consecutive map times were plotted against chan- 
ges in [ol for the same period on a scattergram. The correlation was 
insignificant. This supported the findings of Schardt (12); however, he 
suggested that a correlation might exist between negative values of FG 
(i.e. frontolysis) and decreases in | Vol . Separate scattergrams were 
plotted using only negative FG values but, again, the result was incon- 
clusive. 

7. Conclusions 

In mountainous and sparse data regions and above the 700-mb level in 
the North American area, the time-space continuity of numerical fronts 
was found to be poor. The conclusions stated below refer to the lower 
troposphere east of the Rocky Mountains: 

(a) Numerical baroclinic zones slope toward the cold air or are 
near ly vertical with a few exceptions. Frontal Q values increase with 
height. ° 

(b) Slopes of numerical fronts compare favorably with those of the 
ESSA analyses in some cases, but differ considerably in others. The 
largest deviations are in areas of weak gradient or sparse data. 

(c) GG9 centers have little vertical consistency. Associations 
are best between 850 and 700 mbs. Centers associated with intense baro- 
clinic zones display continuity in time. The more intense the center, 
the more likely it is to persist. Speed of movement was observed to be 


irregular; the 850 and 700-mb levels displayed the best continuity in 


this respect. 
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(d) Average movement of numerical fronts over a sequence of several 
map times is comparable to ESSA frontal movements, but large deviations 
were observed over a given 12-hour interval. 

(e) Of the levels considered, continuity of frontal movement is 
maintained best at the 700-mb level over continental North America. 

(f) FNWF frontal positions, except when the 10/7 model is used, are 
displaced toward the warm air side of the ''true'! hyperbaroclinic zone. 
This appears due to grid size and method of computing derivatives. 

(g) The data processing method employed by FNWF tends to smooth 
open waves in the thermal field and contributes to the present inability 
of the numerical frontal scheme to display proper ly the occlusion process. 
The CAO analyses display open waves in a manner more closely resembling 
the conventional product. 

(h) Frontal analysis at mandatory levels may be improved in accur- 
acy by displacing the GG9 ridge toward the cold air one-half the dis- 
tance from the indicated position to the +0.05°C/(100 km)? isoline. 

8. Recommendations 

The problems involved with lack of data can be solved only with in- 
creased observations and are not germaine to this discussion, There are 
certain aspects of the numerical analysis problem, however, which may be 
susceptible to improvement under present data conditions. The following 
items are listed as suggested areas of further research: 

(a) Reduce the mesh length, thereby reducing the artificial broad- 
ening of the hyperbaroclinic zones and the displacement of the numerical 


front toward the warm air. 
(b) Investigate the feasibility of using observed dry-bulb temper - 
ature as the basic temperature field. This may aid in retaining thermal 


perturbations well into the occluded stage. 
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(c) Where frontal bands can be identified from satellite data, this 
information should be made available to the computer by injecting ''dummy'' 


data in the form of 9, IVa) or GGO. 
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